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NOMENCLATURE 


"1,2.3 

h l. 2.3 


n 

r 

R c 

t 

u 

u 1,2,3 
v 

w 

X 

y 

z 

e 

V 

£ 1,2,3 
CO 

“1,2,3 


unit vectors in transformed space 
vector length in transformed space 
reference length 
nondimensional time index 
radius 

Reynolds number, Uoo 1/v 
time 

velocity component in x-direction 
transformed velocity component 
velocity component in y-direction 
velocity component in z-direction 
Cartesian coordinate 
Cartesian coordinate 
Cartesian coordinate 
angular displacement 
kinematic viscosity 
coordinate in transformed space 
vorticity 

vorticity components in transformed space 


Subscripts 

( )oo 


value at upper plate 



Superscripts 


( )’ dimensional quantity 

0 vector quantity 

( )* quantity at first fractional step 

( )** quantity at second fractional step 

( ) n quantity at nth time step 

( ) n+1 quantity at (n + I )th time step 
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CALCULATION OF THREE-DIMENSIONAL (3-D) INTERNAL FLOW BY 
MEANS OF THE VELOCITY-VORTICITY FORMULATION ON A 

STAGGERED GRID 


Paul M. Stremel 
Ames Research Center 


SUMMARY 


A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) 
enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the 
calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to 
accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the 
non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled 
implicit technique. The solution is calculated on a body fitted computational mesh incorporating a 
staggered grid methodology. In the staggered grid method, the three components of vorticity are 
defined at the centers of the computational cell sides, while the velocity components are defined as 
normal vectors at the centers of the computational cell faces. The staggered grid orientation provides 
for the accurate definition of the vorticity components at the vorticity locations, the divergence of 
vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is 
obtained by utilizing a fractional step solution technique in the three coordinate directions. The 
boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. 

The method provides for the non-iterative solution of the flow field and satisfies the conservation ot 
mass and divergence of vorticity to machine zero at each time step. To test the method, the calcula- 
tion of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow 
in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability 
of the method to predict the flow in arbitrary cavities, results will be shown for both cubic and 
curved cavities. 


INTRODUCTION 


The numerical prediction of vortex-dominated flow is paramount to the understanding ot the 
flow about aircraft configurations. This is especially important to the design analysis of rotorcratt 
when considering download, the force exerted on the vehicle due to the rotor-wake interaction with 
the fuselage components. Download limits helicopter performance in hover and is a significant 
problem in the design of tilt-rotor configurations, where the lifting wing is immersed in the rotor 
wake. Because the download caused by the rotor wake severely limits the hover performance of tilt- 
rotor configurations, a method for accurately predicting tilt-rotor download would provide for the 
design of configurations with improved hover performance. 



A two-dimensional (2-D) method has previously been developed to calculate the flow about bluff 
bodies, reference 1. Results have also been obtained for airfoils with and without a deflected flap at 
-90 deg incidence, reference 2. Additionally, the effect of Reynolds number and turbulence have 
been computed for the XV- 15 wing airfoil with and without a deflected flap, reference 3. The results 
of reference 3 indicate that the flow field solution is highly Reynolds number and turbulence depen- 
dent. Excellent correlation between prediction and test were obtained when matching the test 
Reynolds number and incorporating the Baldwin/Barth turbulence model, reference 4. This correla- 
tion provides confidence in using the current method as a tool to further investigate the reduction of 
drag on airfoils at -90 deg incidence. 

The validated 2-D computational method has also been applied to calculate the influence of 
upper- or lower- surface fences on airfoil aerodynamics, reference 5. In particular, the flow about an 
XV- 15 airfoil with a 30 percent trailing edge flap deflected 60 deg at -90 deg incidence was con- 
sidered. The flow is calculated for a Reynolds number of one million while modeling turbulent flow. 
The results of that investigation indicate that significant reductions in drag are obtained with the 
inclusion of fences. In particular, a 35 percent drag reduction, with respect to the basic airfoil value, 
was achieved for a lower-surface fence located at the airfoil leading edge. 

The ability of the 2-D method to accurately compute the flow about a complex geometry normal 
to the free-stream flow and the direct extension of the method to 3-D analysis are the basis for this 
paper. In order to predict download, a method is required that can, not only, compute the flow about 
complex 3-D bodies, but also can accurately predict the wing-base pressure. The 3-D extension of 
the 2-D analysis promises to be such a method. 

Prior to solving the 3-D external flow problem a simpler problem is addressed to validate the 
governing equations and the solution technique. The flow in a driven cavity is considered to test the 
method. The development of the method is presented in the next section followed by the application 
of the method to the driven cavity problem. 


PROBLEM FORMULATION 


The flow field is modeled by the velocity/vorticity form of the unsteady, incompressible Navier— 
Stokes equations. The nondimensional governing equations in Cartesian coordinates are written for 
the continuity equation. 


V • u = 0 


0 ) 


and for the vorticity transport equation. 


co t +V x(cbxu) = V 2 co/ R e (2) 

with V 2 = d 2 ( )/d x 2 + 3 2 ( )/dy 2 + d 2 ( )/3z 2 , where (x,y,z) are the Cartesian coordinates, R e is the 
Reynolds number, and t is the time. The vorticity, 0), is defined by 


-> 



a) = V x u 


(3) 


The nondimensional variables are written 

x = x' / l,y = y’ / l,z = z' / / 
m = w ' / U^v = v' / U^w = vv' / 

(0 = co'/{U BO /l),t = t'/{l/U oa ),R e =UJ/v 


where 


u, v, and w are the Cartesian components of the velocity, and 

1 = reference length 

R c = Reynolds number based on I 

Uoo = free-stream velocity 

v = kinematic viscosity 

When the Cartesian equations are transformed into generalized-orthogonal-curvalinear coordi- 
nates, the governing equations become as follows (see ref. 6 for details). 


For the continuity equation. 


V • (7 = 




d(h 2 hjUi) d(h,h 3 u 2 ) d(h,h 2 u 3 ) 

Hi + H 2 + Hi 


= 0 
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and for the definition of vorticity. 
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where (^i,^2^3) are the transformed coordinates, u ( , u 2 , and u 3 are the orthogonal velocity compo- 
nents in the transformed coordinates, and hj, h 2 , and h 3 are the vector lengths. 
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The vorticity transport equation, equation 2, can be rewritten as: 

<y ; +M»V6>-u>»Vii = / R e 


(6) 


This form offers some unique issues regarding the dot product on the left-hand side and the 
Laplacian operator on the right-hand side of the equation. Because both the gradient and the 
Laplacian operator act on a vector quantity, the dependence of the vector components and the unit 
vectors must be considered. The expression for 5 • Vai becomes: 
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Where, V( ) = 


!<?()_ 1 <?( ), _ 1 d( ) _ 

— a, + — flb + — 

A, h 2 Hi ‘ Aj 


. Details can be found in Appendix 1 . 


The expression for the Laplacian of the vorticity vector on the right-hand side of equation 6 is 
extremely long and can also be found in Appendix 1 . 

The boundary conditions for the transformed governing equations at the enclosure surface are 
calculated from the no-slip condition as 


M| = «2 = M3 = 0 (7) 

except at the upper surface of the enclosure at which the velocity is specified. The surface vorticity is 
calculated from equation 5. 
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The flow is started impulsively. Therefore, at t = 0, the velocity on the upper surface of the 
enclosure is set to the prescribed boundary condition and the velocity everywhere inside the enclo- 
sure is set equal to zero. The vorticity on the upper surface of the enclosure is calculated from 
equation 5 and is set equal to zero everywhere else in the enclosure. The surface vorticity is calcu- 
lated implicitly as part of the solution after the impulsive start. The computation is advanced until a 
desired time has been reached or until the flow has demonstrated stable periodic flow within the 
enclosure. 


NUMERICAL METHOD 


The solution is obtained by solving the finite-difference representations of the governing equa- 
tions on a computational mesh. The grid is body-fitted to the interior of the 3-D enclosure. 

In the staggered-grid method, the flow-field variables are not defined at the mesh nodes only. 
Rather, the components of vorticity are defined at the mid-points of the mesh cell sides, and the 
orthogonal flow-field velocity components are defined at the centers of the mesh cell faces. The vor- 
ticity and velocity components on the staggered grid are depicted in figure 1 . The staggered-grid ori- 
entation of the variables provides for the conservation of vorticity at the mesh nodes and the solution 
of the continuity equation at the mesh cell centroids. The vorticity and flow-field velocities are cal- 
culated by a fully coupled implicit technique on the staggered mesh. The coupled method solves for 
the vorticity and velocity components by means of a block-tridiagonal inversion for fractional steps. 
A representation of the fractional step method is presented in Appendix 2. Each fractional step repre- 
sents a computational sweep in one of the coordinate directions. These computational sweeps are 
depicted in figures 2^1. 

In figure 2, the computational sweep for coordinate is shown. The governing equations for the 
vorticity and velocity components are selected to take into account the spatial derivatives in 2,j . With 
this in mind, the conservation of vorticity is solved for CO] . The vorticity components (1)2 and 0)3 are 
solved using the second and third components of the vorticity transport equation, equation 6. The 
continuity equation, equation 4, is solved for U|, and the velocity components U2 and U3 are solved 
from the third and second components of vorticity, respectively, in equation 5. This allows for the 2,i 
derivatives of U2 and U3 to appear in the governing equation for each velocity component. Then 
for U2, 


and for U3, 
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In figure 3 , the computational sweep for coordinate ^2 is shown. The governing equations for the 
vorticity and velocity components are selected to take into account the spatial derivatives in ^2- Now, 
the conservation of vorticity is solved for 0)2. The vorticity components a>i and 0)3 are solved using 
the first and third components of the vorticity transport equation, equation 6. The continuity equa- 
tion, equation 4 , is solved for U2, and the velocity components ui and U3 are solved from the third 
and first components of vorticity, respectively, in equation 5 . This allows for the ^2 derivatives of U| 
and U3 to appear in the governing equation for each velocity component. Then for U|, 


and for U3, 


ft >3 = 


h A 


^(Vg] 

Hi 
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Hi J 
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In figure 4 , the computational sweep for coordinate ^3 is shown. The governing equations for the 
vorticity and velocity components are selected to take into account the spatial derivatives in £3. Now, 
the conservation of vorticity is solved for (1)3. The vorticity components coj and 0)2 are solved using 
the first and second components of the vorticity transport equation, equation 6. The continuity equa- 
tion, equation 4 , is solved for U3, and the velocity components uj and U2 are solved from the second 
and first components of vorticity, respectively, in equation 5 . This allows for the £2 derivatives of U] 
and U3 to appear in the governing equation for each velocity component. Then for ui. 
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and for U2, 
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Hi 
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The boundary conditions are calculated implicitly as part of the solution. After the inversion of 
these block-tridiagonal systems, the solution for the present time-step is obtained. The solution is not 
iterated and the continuity equation is satisfied to machine zero. The solution is then advanced by 
updating the flow-field variables and solving the flow for the next time-step. The solution is termi- 
nated at a desired time or after sufficient time has elapsed for demonstration of the flow-field 
periodicity. 
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RESULTS 


In order to test the method, the solution to the driven cavity problem was computed. All solutions 
were computed for a Reynolds number of 100. This geometry was considered to test the methodol- 
ogy and not to analyze the flow of driven cavities. The current formulation provides for the solution 
of 3-D, incompressible flow for arbitrary geometries and not constrained to Cartesian coordinates. 
The driven cavity problem is applied for a cubic cavity with constant coordinate spacing and with a 
cosine spacing, which clusters the grid lines at the boundaries. The results are compared to those 
computed in reference 8 for a cubic cavity of constant grid spacing in Cartesian coordinates. A 
curved cavity is .also considered, again with a constant coordinate spacing and a cosine spacing for 
the grid coordinates. The cross-flow geometry for the four grids considered are shown in figures 5-8. 
The 3-D grids are constructed by stacking the cross-flow geometry. The velocity of the upper surface 
of the cavity is specified and held constant. The solution is started impulsively from r.est and is con- 
tinued until a steady solution is reached. 


Cubic Cavity - Constant Grid Spacing 

The driven cavity flow for the geometry of figure 5 is computed. For this solution, x is synony- 
mous with 4i, y with ^2- and z with £3. The upper plate, the x-y plane at maximum z, is held at a 
constant velocity, Uo° = 1 -0. The vorticity contours at the mid-plane locations are shown in figure 9. 
The vorticity, shown in figure 9(a), demonstrates the symmetric nature of the flow field with 
opposite vortices at the center of the plane and flow up the walls at y = 0 and y = 1.0. Additionally, 
secondary vortices form on the bottom of the cavity and on the upper plate. The ^2 vorticity, shown 
in figure 9(b), develops similar to the vorticity for a 2-D solution of the driven cavity problem. The 
£3 vorticity, shown in figure 9(c), demonstrates the symmetric nature of the cavity flow with sec- 
ondary vortices on the upstream wall. The vorticity contours of figure 9 are similar to the vorticity 
contours presented in reference 8. 


Cubic Cavity - Cosine Grid Spacing 

The driven cavity flow for the geometry of figure 6 is computed. This is similar to the previous 
solution. However with the increased grid resolution at the wall, additional flow features in the form 
of secondary vortices are resolved. Here again, x is synonymous with £1. y with ^2- and z with ^3 and 
the upper plate is held at a constant velocity, Uoo = 1 .0. The vorticity contours at the mid-plane loca- 
tions are shown in figure 10. The ^1 vorticity, shown in figure 10(a), demonstrates the symmetric 
nature of the flow field with opposite vortices at the center of the plane and flow up the walls at 
y = 0 and y = 1 .0 and the enhanced flow features, such as, secondary vortices on the bottom of the 
cavity and on the upper plate. The t,2 vorticity, shown in figure 10(b), develops similar to the 
vorticity in the 2-D solution of the driven cavity problem and the previous vorticity shown in 
figure 9(b). Vortices in the lower corners of the cavity are clearly predicted. The ^3 vorticity, shown 
in figure 10(c), demonstrates the symmetric nature of the cavity flow with secondary vortices on the 
upstream and downstream wall. The vorticity contours in figure 10 are similar to those of figure 9 
and the vorticity contours presented in reference 8. 
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Curved Cavity - Constant Grid Spacing 


The driven cavity flow for the geometry of figure 7 is computed. For this solution, r, or radius, is 
synonymous with £j, y with £2. and 0 with ^3. The upper plate, the x-y plane at maximum z, is held 
at a constant velocity, Uoo = 1 .0, see figure 7. The vorticity contours at the mid-plane locations are 
shown in figure 1 1. The vorticity, shown in figure 11(a), again demonstrates the symmetric nature 
of the flow field with opposite vortices at the center of the plane and flow up the walls at y = 0 and 
y = 1.0. Secondary vortices shown for the cubic cavity, figure 9(a), are not present here owing to the 
curvature and increased height of the cavity, both of which tend to stagnate the flow at the bottom of 
the cavity. The ^2 vorticity, shown in figure 1 1(b), develops similarly to the vorticity shown in. 
figure 9(b) for the cubic cavity. However, the vorticity develops further down the downstream wall 
due to the curvature of the cavity. The £3 vorticity, shown in figure 1 1(c), demonstrates the symmet- 
ric nature of the cavity flow with secondary vortices on the upstream wall. 


Curved Cavity - Cosine Grid Spacing 

The driven cavity flow for the geometry of figure 8 is computed. Here r, or radius, is synony- 
mous with £1, y with £2* and 0 with ^3. The upper plate, the x-y plane at maximum z, is held at a 
constant velocity, Uoo = 1.0. The vorticity contours at the mid-plane locations are shown in figure 12. 
The vorticity, shown in figure 12(a), again demonstrates the symmetric nature of the flow field 
with opposite vortices at the center of the plane and flow up the walls at y = 0 and y = 1.0. The sec- 
ondary vortex structure shown on the upper plate is better defined, when compared to that of the 
constant grid spacing, figure 1 1(a). The ^2 vorticity, shown in figure 12(b), develops similar to the 
vorticity shown in figure 1 1(b). However, comer vortices can be seen in the bottom of the cavity and 
increased rollup of the primary vortex in the center of the cavity. The £3 vorticity, shown in 
figure 12(c), demonstrates the symmetric nature of the cavity flow with secondary vortices on the 
upstream and downstream walls. 


CONCLUSIONS 


A method has been developed to calculate accurately the viscous flow in 3-D enclosures. The 
method provides for the non-iterative solution of the incompressible Navier-Stokes equations by 
means of a fully coupled implicit technique. To demonstrate the method, the calculation of simple 
driven cavity flows has been considered. The driven cavity flow is defined as the flow in an enclo- 
sure driven by a moving upper plate at the top of the enclosure. The intent was not to present a study 
of cavity flows, but to determine the ability of the method to predict arbitrary internal 3-D flow. 
Therefore, results were shown for both cubic and curved cavities with constant and varying mesh 
spacing. 

The predicted vorticity contours for a cubic cavity with constant mesh spacing were in close 
agreement with the contours shown in reference 8. The method predicted the emergence of the pri- 
mary vortices, as well as, secondary vortices. Similar results were shown for a cubic cavity with a 
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cosine spacing for the grid points. The cosine spacing of the grid points concentrates the mesh points 
near the boundaries. With the cosine spacing, the primary and secondary vortices were enhanced. 

The predicted vorticity contours for a curved cavity with constant mesh spacing were similar to 
those for the cubic cavity, and were again enhanced when the solution was computed on a curved 
cavity with a cosine distribution for the grid points. 

The comparison of the present results with other computations for the cubic cavity flow and the 
comparison of the cubic cavity flow with the curved cavity flow provide credibility to the methodol- 
ogy. Considerable work is required to extend the method to external flows, to incorporate turbulence 
models and to apply the method to realistic problems. 


9 




APPENDIX 1 


This Appendix presents the development of the dot-product, divergence, and LaPlacian operator of a 
vector quantity and are integral to the development of the generalized governing equations in three- 
dimensional (3-D) space. 

The expressions in this Appendix are written for a general vector quantity. 


F = F|«| + +F 2 a 2 + F 2 a 2 


The dot product of a vector with the gradient of another vector is written, 

ii • VF = n • V(F|rt| + F\fi 2 + F^) 


or, 

n • VF = n • VF|«i + n •VF 2 a 2 + w • V F 2 a 2 + F|(/7 • V«|) + F 2 (/7 • V« 2 ) + * V 03 ) 


The general expressions for the dot product and the gradient are, respectively. 
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(he dot product of the unit vectors become. 
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Then the dot product becomes, 


or. 
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The Laplacian of a vector is as follows. 
The Laplacian of a scalar is written. 
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For which the Laplacian for a vector quantity is written. 
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After evaluation of the above terms using the product rule, incorporating the derivatives of the unit 
vectors presented above and considerable rearranging, the Laplacian of a vector quantity is written, 
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h\h 3 V<9^| J h x h 2 \d^ x ) h x h 2 \d^ 2 J h x h 3 ^<9£ 3 J 

+ 1 ^ h 3 dh\ dF 2 | ^ d ( Ih, dh } ' _. 

h\h 2 h 3 /i[ <9 £ 2 <9<^i ^2 J_ 

, 1 2 Ih ^1 ffijt | /- - 

h\h 2 h 3 ~ h\ <9 £ 3 <9^| <9^| [ /i| <9<^3 

I /h dh, «9F, r9 f h 2 dho ) - 

_j_ — 2 — : — — jT 2 — : ~ £/| 

/ij/l 2 /l 3 ^2 ^^2 *^>2x^2 ^>1 / 

1 \ ( dh 3 <?/l 3 ^ 1 ^ dlh dll2 > _ 

h\h 2 h 3 h 3 v <9Ci <9^| > h 2 y d £ 3 <9^| 

1 ^ h 2 dh 3 dF 3 d f h 2 dh 3 ' _ 

h\h 2 h 3 ~ ~ h 3 dt ; | <9^3 9^3 ^3 <9£, ^ 
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d f h 2 h 3 c)F 2 ] d ( /ij /»3 dF 2 1 d f h ] h 2 dF 2 

h\h 2 h 2 /, ( <9/13 Y h 2 ( dh 2 \ h\ c f dh 2 1 h 2 r f dh, 1 

^2^3 1^2/ h^2 ) hohy ) h\ho v<9^2 ) 

, ' r,», ^ <#■, r <? f *» ^ ii, 

*1*2*1 ft, d^2 1 *^l 1*1 ^22 J 

+ 1 2 ^3 dh 2 dF\ + p d h 2 dh 2 ' _ 

/i| h 2 h 2 h 2 <9£, r^ 2 1 <9^2 v ^2 ^Si > 2 

1 /?> dh 2 dF-, d f lh dh 2 1 _ 

+ 2— 2. — i +F 2 — f- ci 2 

h^hoh^ ho <9^3 <9^2 <9<^2 y hn <9^3 -j 

1 1 f dh 2 dh 2 ' 1 ^ dh^ dh\ ' 

h\h 2 h 2 h 2 y o^i <9^2 > v <9s3 ^2>. 

| i f —(^ ^3 1 - 

h\h 2 h 2 h 2 d% 2 d^ 2 v /*3 <9<^ 2 , 


’ d 

hoh-\ 

<9F3 

<9 

-f - • - - 

/; ( /<3 r9F 3 ' 

d 

+ 

/j,/l2 <9^ ^ 


l 


J ^2 

v hi d £ 2 / 

) ^3 

v /«3 <9^3 j 


W 3 h 2 „ f dh-\ V /j, „ f <9/^ V //, r f dh 2 f /» 2 r f d/i, f 

/)|/l 3 /l 2 /?3 d^ 2 J ^2^3 ^3 J ^1^3 <9<53 > 

_1 ^ ^<9/^ _y_f ^2_^]_1 - 

h\h 2 h 2 *” ^3 <?£| <9<^i ^ /i] <9& > _ 

I i f dfh dli-, 1 i f <9/;, <9/i| 1 _ 

— /r — _ — - + — F 2 — 1 L «3 

h\h 2 h 2 h 2 dc i \ ^3 J /?| <9(^2 / 

I /i, dh, dF-, f <9 f /i, dlh ^ _ 

+ 2— 1 -- F 2 — 1 «3 

h\ll 2 h 2 h 2 <9(^3 i9i ^ 2 <9(^2 h 2 <9^3 /J 
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| 1 2 ^1 ^3 I f f h\ dh 3 ^ _ 

*1*2*3 L *3 ^3 2 ^,U^Tjj'’ 
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APPENDIX 2 


This Appendix presents the development of the fractional-step method of reference 7 for three- 
dimensions. Consider the Cartesian form of the vorticity transport equation in three dimensions, 
without considering the vortex-stretching terms. The equation is written, 

dtO d{uto) d(vto) d{w(t)) _ 1 d to + d to + d CO 
dt dx dy dz Re\dx 2 dy 2 dz 2 > 

The fractional step method is written , for the three fractional steps, as 
for the x-coordinate step, 

a)*-(o n 1 d(uto)* | d{u(o) n ^ | d(v(o) n | <9(wtu) n _ J_ l_ d^to_ | d^to^ | d^to^ | d^to^_ 

At 2 ^ dx dx j dy dz Re K 2 y ^x 2 dx 2 } dy 2 dz 2 ) 

for the y-coordinate step, 

to** - to 11 1 d(uto) + d{uto) n + I d(vco) + d(vto) n | d{w(6) n _ 

At 2 K dx dx j 2 v dy dy } dz 

1 l( d 2 tO* d 2 tO n ) 1 ( d 2 (0** d^o/^ d^to^ 

^[2[d?~ + dx 2 j^{ d? dy 2 ) dz 2 , 

and for the z-coordinate step, 

+ l( d{u(o)* 7 ( djvtof | l( d{wto) n + 1 1 d(w»') n \._ 

At 2 ^ dx dx j 2 y dy dy J 2[ dz dz. } 

i(l(d 2 tO* d 2 (O n ) lid 2 CO** d 2 0) n 1 l( d^0)^^_ d^O^_ V 
at 2 ) + 2( d? + ~d?~)^{ dz 2 + dz 2 J, 

Then 


I * n\ 
to -to 

2 U 

At 


d(uto) + d(uto )' 1 
dx dx > 


d{vto) n + ^d{ wto) n 
dy dz 


. ( ( -,2 * -.2 n\ 

I d to d to 

Ro y y dx 2 dx 2 ) 


cTto 

dy 2 


+ 2 



(a>** -ffl*) 
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yields, 


yields, 


, ( d(vco)** <9(vw)'0 l(d 2 aT d 2 co n ^ 


At 


-+ 


dy 


d\’ 


Re 


dv‘ 


dy‘ 


(O n+] - CO** 


(co” + , -co**) 


At 


+ 


<?(wft))” +/ djwwf ^ 
dz dz 


Re 


d 2 co n+ 1 d 2 co n 1 


dz‘ 


dz‘ 
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Figure 1. Schematic for staggered-grid formulation. 
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Figure 5. Geometry for cubic cavity, constant spacing, 17 x 17 x 17. 
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Figure 6. Geometry for eubic cavity, cosine spacing, 17x17x17. 
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Figure 7. Geometry for curved cavity, constant spacing, 17 x 17 x 17. 
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Figure 8. Geometry for curved cavity, cosine spacing, 17 x 17 x 17. 
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